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We derive here a linear elastic stochastic description for slow crack growth in heterogeneous 
materials. This approach succeeds in reproducing quantitatively the intermittent crackling dynamics 
observed recently during the slow propagation of a crack along a weak heterogeneous plane of a 
transparent Plexiglas block [Mal0y et al., PRL 96 045501]. In this description, the quasi-static 
failure of heterogeneous media appears as a self-organized critical phase transition. As such, it 
exhibits universal and to some extent predictable scaling laws, analogue to that of other systems 
like for example magnetization noise in ferromagnets. 

PACS numbers: 62.20.mt, 46.50. +a, 68.35. Ct 



Driven by both technological needs and the challenges 
of unresolved questions in fundamental physics, the effect 
of materials heterogeneities onto their failure properties 
has been extensively studied in the recent past (see [l| 
for a recent review). So far, many efforts have been fo- 
cusing on the morphology of fracture surfaces (see [3| for 
a review). In particular, crack surface roughness was re- 
cently shown to exhibit anisotropic morphological scaling 
features 0, 0] that could be understood for brittle mate- 
rials [3]. 

Here, we will focus our study on the dynamics of 
cracks. In heterogeneous materials under slow external 
loading, this propagation displays a jerky dynamics with 
seemingly random sudden jumps spanning over a broad 
range of length-scales [3] . Such a complex dynamics also 
called "Crackling Noise" has been also suggested from 
the acoustic emission accompanying the failure of various 
materials Q and the seismic behaviour accompanying 
earthquakes Q characterized by power-law energy distri- 
bution with an exponent around 1.3 1.5. These distri- 
butions can be qualitatively reproduced in phenomeno- 
logical models like Fiber Bundle Models (FBM) (j| (resp. 
Random Fuse Models (RFM) [l(J) which schematizes ma- 
terials as a set of brittle fibers loaded in parallel with ran- 
dom failure thresholds and a rule for load redistribution 
after each failure event (resp. as a network of electrical 
fuse with constant resistance and randomly distributed 
threshold). However, these simple approaches yield to an 

exponent around 1.9 2.5 significantly higher than the 

experimental observations. Moreover, they rely on im- 
portant simplifications which makes quantitative com- 
parisons with experiments difficult. 

We will demonstrate here that this crackling dynam- 
ics can be fully reproduced through a stochastic descrip- 
tion rigorously derived from Linear Elastic Fracture Me- 



chanics (LEFM) [l2j extended to disordered materials. 
In particular, this model succeeds in reproducing quanti- 
tatively the intermittent crackling dynamics observed re- 
cently during the steady slow crack growth along a weak 
heterogeneous plane within a transparent Plexiglas block 
1 11 ] - In this description, quasi-static failure of heteroge- 
neous brittle elastic media can be interpreted as a self- 
organized critical dynamic phase transition and - as such 
- exhibits universal behaviors. We will then show how one 
can use Universality and previous calculations performed 
on different systems belonging to the same universality 
class, here the motion of domain walls in disordered fer- 
romagnets, to derive predictive laws for the failure of ma- 
terials. 

Theoretical description: LEFM is based on the fact 
that - in an elastic medium under tensile loading - the 
mechanical energy G released as fracture occurs is en- 
tirely dissipated within a small zone at the crack tip. 
Defining the fracture energy T as the energy needed to 
create two crack surfaces of unit area, under quasi-static 
condition, we assume that the local crack velocity v is 
proportional to the excess energy locally released: 



-v = G-T 



(1) 



where fi is the effective mobility of the crack front. 

In a homogeneous medium, r = Tq is constant and 
an initially straight crack front will be translated with- 
out being deformed. LEFM allows the determination of 
the energy released G in any geometry. In particular, for 
the experimental configuration chosen in [ll| , considered 
here and depicted in Fig.l where a crack grows at the in- 
terface between two elastic plates by lowering the bottom 
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FIG. 1: Crack growth along the inhomogeneous interface be- 
tween two elastic plates loaded according to the geometry 
used in 11 ill: Sketch and notations. 



plate at constant velocity V, G is given by [131 ] : 

i 3 (A +vtr 



G°(f(t)) = -E6 
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(2) 



where c + f(t) is the instantaneous crack length, E the 
Young modulus, S the thickness of the lowered plate, Ao 
the initial opening and Co the initial crack length. Then, 
considering a slow driving velocity such as Vt <C Ao, 
one can show from Eqs. II 121 that after a short transient 
regime, the crack front propagates at a constant velocity: 



f(t) — v m t with v n 
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In a heterogeneous material, defects induce fluctua- 
tions in the local toughness: T(x,y) = T°(l + ri{x,y)) 
where the noise term r)(x,y) captures these fluctuations, 
and x and y denotes the coordinates in the propagation 
and crack front directions respectively. These fluctuations 
induce local distortions of the crack front f(x, t), which 
in turn generate local perturbations in G [l4j]. To linear 
order in /, one can show that [li| : 

G(x,f(x,t))v G°((f(x,t)) x ) 



+iG»((/(M)Un&^)rf,' 

(4) 

And finally, by replacing this expression in Eq. [TJ using 
the expression of G° for the homogeneous case, and by 
introducing the quantities defined in Eq. [3J we get: 

FH, {/}) + h C fiX '$:^ dx> - V (x, f(x, t)) 
where F(t,{f}) = ±(v m t- < f(x,t) > x ) 

(5) 

Strictly speaking, this equation describes the interfacial 
crack growth according to the geometry depicted in Fig. 
1. But the very same equation - with different prefactors 
in the expression of F(t, {/}) = ^ (v m t- < f(x,t) > x ) 
- would have been obtained for the quasi-static stable 
crack growth of the in-plane component of the crack 
front within a three-dimensional solid independently of 
the tensile loading conditions and the system geometry. 



Variants of this equation with constant force F have 
been extensively studied in the past to model crack prop- 
agation in solids but also to describe other systems 
as diverse as interfaces in disordered magnets [3, Qj| 
or contact lines of liquid menisci on rough substrates 



2d l2ll |. In this case, the interface remains stationary, 



pinned by the heterogeneities, unless this constant force 
F is larger than a threshold value F c . A key feature of 
these systems is that the so-called depinning transition at 
F c belongs to the realm of collective critical phenomena 
characterized by universal scaling laws j^. In particular, 
at F c , the interface moves through scale free avalanches, 
both in space and time. 

Equation [5] denotes a rather different situation where 
the effective force F(t, {/}) is not constant anymore, but 
given by the difference between the mean front position 
and the one that would have been observed within the 
homogeneous case: When F(t, {/}) < F c , the front re- 
mains pinned and F(t, {/}) increases with time. As soon 
as F(t,{f}) > F c , the front propagates, F(t, {/}) in- 
creases, and, as a consequence, F(t, {/}) is reduced until 
the front is pinned again. This retroaction process keeps 
the crack growth close to the depinning transition at each 
time and, within the limit v m — > and Co — > oo the sys- 
tem remains at the critical point during the whole prop- 
agation, as for self-organized-critical systems [23]. 

Spatio-temporal intermittent dynamics - Predictions of 
this stochastic description are now confronted to the ex- 
perimental observations reported in [ll|. Using a fourth 
order Runge-Kutta scheme, Eq. [5] is solved for a front 
f(x,t) propagating in a 1024 x 1024 uncorrelated 2D 
random gaussian map rj(x,y) with zero average and unit 
variance [Note that the crackling dynamics statistics de- 
tailed thereafter do not depend on the precise distribu- 
tion of r] as far as it remains short-range correlated] . The 
parameter ^lTq was set to unity while the two parame- 
ters co and v m were varied from 2.5 to 250, and 10 -2 to 
5 x 10" 1 , respectively. In order to characterize the scal- 
ing features of the crack front local dynamics, we adopt 
the analysis procedure recently proposed in [ll| and com- 
pute at each point (x, y) of the recorded region the time 
w(x, y) spent by the crack front within a small 1 x lpixel 2 
region as it passes through this position. A typical gray- 
scale image of this so-called waiting time map are pre- 
sented in Fig. 2a. The numerous and various regions of 
gray levels reflect the intermittent dynamics, and look 
very similar to those observed experimentally (Fig. lb 
of (ill). From the inverse value of waiting-time map, we 
compute a spatial map of the local normal speed veloc- 
ity v(x, y) = l/w(x,y) of the front as it passes through 
the location (x,y). Avalanches are defined as clusters of 
velocities v larger than a given threshold v c = Cv where 
v denotes the crack velocity averaged over both time and 
space within the steady regime. Their area is power-law 
distributed with an exponent tq = — 1.65 ± 0.05 (Fig. 
2c) . These clusters exhibit morphological scaling features 
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since their width, - measured along the direction of crack 
growth is shown to go as a power-law with respect to 
their length - measured along the direction of crack front 
- with an exponent H = 0.65±0.05 (Fig. 2d). All these re- 
sults are in perfect agreement with those measured exper- 
imentally in [lH. Contrary to what conjectured in 11 1, 
H is significantly different from the value of the rough- 
ness exponent ( ~ 0.39 expected for the interface at the 
depinning point of Eq. [5] [3, E3| ■ It is worth to mention 
that recent experimental results [17J report a roughness 
exponent C — 0.35 at large scale, in agreement with this 
theoretical prediction. 

To complete the characterization of the avalanche 
statistics, we measure a new observable, the avalanche 
duration, - defined as the difference between the times 
when the crack front leaves and arrives to the consid- 
ered avalanche cluster (Fig. 2b). The distribution of the 
avalanche duration is plotted on Fig. 2e. For durations D 
smaller than the average avalanche duration < D >, this 
distribution clearly depends on v m , cq and the clipping 
threshold C while for D 3>< D > all the distributions 
seem to collapse onto a single power-law behaviour char- 
acterized by an exponent a> = —2.1 ± 0.2 (Fig. 2e). 
Finally, the mean avalanche duration is found to go as 
a power law with the mean avalanche size, characterized 
by an exponent 7 = 0.4 ± 0.05 (Fig. 2f). 

Spatially- averaged dynamics - As for other critical sys- 
tems, very different physical systems that belong to the 
same universality class will display similar scaling be- 
haviours with same scaling exponents. In particular, it 
was shown that the Barkhausen noise fljl . [l9| accompa- 
nying the motion of domain wall driven by external mag- 
netic field through soft disordered ferromagnets can be 
described by an equation with shape similar to Eq. [5j In 
this respect, we analyse the global crack front dynamics 
by computing the spatially averaged instantaneous ve- 
locity?;^) =< ^(x, t) > x of the crack interface, as in 
13, El- We observe once again a jerky dynamics really 
similar to Barkhausen noise [3, 0|- We then impose a 
given reference level v c = Cmax(i)) and define bursts as 
zones where v(t) is above v c . The duration T of a given 
burst is defined as the interval within two successive in- 
tersections of v(t) with v c while the size S is defined 
as the integral of v(t) between the same points. For the 
Barkhausen noise, this analysis leads to power-law dis- 
tributions P(T) cx T- a and P(S) cx S~ T [H, Hi, and 
a power-law relation between the duration and size of 
the avalanches T cx S a , with critical exponents that can 
be predicted using Functional Renormalisation Group 
(FRG) calculations OHO] leading to r « 1.25, a « 1.43, 
and a m 0.58. We show on Fig. 3b, c and d. that all these 
scaling relations as well as the values of the exponents are 
found to be in good agreement with the ones observed 
here for the average interfacial crack growth. 

Conclusion - We have derived a description for pla- 
nar crack growth in a disordered brittle material which 
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FIG. 2: (a) Typical gray scale map of the waiting time matrix 
w(x,y) obtained from the solution of Eq. [5] with Co = 4000 
and Vm — 0.01. (b) Spatial distribution of clusters correspond- 
ing to velocities v(x,y) > Cv with C = 4. The quantity v 
refers to the crack velocity averaged over both time and space 
within the steady regime. The clusters duration is given by 
the clusters color according to the colorscale given in inset. 
The distribution of the cluster area A, the scaling between 
the width L y and the length L x of clusters, the distribution 
of the cluster duration D, and the scaling between D and 
A are plotted in (c), (d), (e) and (f), respectively. The vari- 
ous symbols corresponds to various values of v m , Co and C as 
specified in the inset of Fig. e. The straight lines correspond 
to P(A) oc A~ T °, L y cc Lf , P(D) oc D~ a > and Da A 7 with 
to = -1.65, H = 0.65 ± 0.05, a> = -2.1 and 7 = 0.4. 



succeeds to capture the statistics of the intermittent 
crackling dynamics recently observed experimentally [11|. 
In particular, we have shown that material failure ap- 
pears as a critical system, where the crack front pro- 
gresses through scale-free avalanches signature of a dy- 
namic depinning phase transition. As for other critical 
systems, microscopic and macroscopic details do not mat- 
ter at large length and time scales and this simple Lin- 
ear Elastic Stochastic description contains all the ingre- 
dients needed to capture the scaling statistical proper- 
ties of more complex failure situations. This conjecture 
is strongly supported by a recent analysis [25j that has 
shown that the laboratory experiment described here of 
crack gowth within Plexiglas under tension share many 
characteristics with geological faults - although this lat- 
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FIG. 3: (a): Typical time evolution of the spatially averaged 
crack front velocity v(t) =< v(t,x) > x . The dash-dot hori- 
zontal line corresponds to v r = Cmax(w) with C = 0.1 used 
to define the bursts. Distribution of the normalized burst du- 
ration T, normalized burst size S, and the scaling of T with S 
are plotted in (b), (c) and (d), respectively. The various sym- 
bols corresponds to various values of v m , Co and C as specified 
in the inset of Fig. b. In these three graphs, the straight lines 
correspond to P(T) oc T~ a , P(S) <x S~ T and T tx S a with the 
critical exponents a = 1.43, r = 1.25, and a = 0.58 predicted 
by FRG approach. 



ter results from external shear loading - since they ex- 
hibit seismicity catalogs with similar satistical features: 
Time of occurrence, epicentre location and energy pa- 
rameter. FRG calculations performed on a different sys- 
tem belonging to the same universality class, - the mo- 
tion of domain walls in disordered ferromagnets have 
been used to predict scaling laws for the statistics of the 
spatially averaged instantaneous crack velocity. Exper- 
imental checks of these predictions are currently under 
progress Let us note finally that this description has 
been derived within the quasi-static approximation. To 
understand how to include the dynamic stress transfers 
through acoustic waves occurring as a dynamically grow- 
ing crack is interacting with the material disorder [26j will 
represent interesting challenges for future investigations. 
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